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Abstract 

Increased precipitation is one projected outcome of climate change that may enhance the discharge of freshwater to 
the coastal zone. The resulting lower salinity, and associated discharge of both nutrients and dissolved organic 
carbon, may influence food web functioning. The scope of this study was to determine the net outcome of increased 
freshwater discharge on the balance between auto- and heterotrophic processes in the coastal zone. By using long- 
term ecological time series data covering 13 years, we show that increased river discharge suppresses phytoplankton 
biomass production and shifts the carbon flow towards microbial heterotrophy. A 76% increase in freshwater 
discharge resulted in a 2.2 times higher ratio of bacterio- to phytoplankton production (P b :P p ). The level of Pb:P p is a 
function of riverine total organic carbon supply to the coastal zone. This is mainly due to the negative effect of fresh- 
water and total organic carbon discharge on phytoplankton growth, despite a concomitant increase in discharge of 
nitrogen and phosphorus. With a time lag of 2 years the bacterial production recovered after an initial decline, further 
synergistically elevating the microbial heterotrophy. Current climate change projections suggesting increased precipi- 
tation may therefore lead to increased microbial heterotrophy, thereby decreasing the transfer efficiency of biomass 
to higher trophic levels. This prognosis would suggest reduced fish production and lower sedimentation rates of phy- 
toplankton, a factor of detriment to benthic fauna. Our findings show that discharge of freshwater and total organic 
carbon significantly contributes to the balance of coastal processes at large spatial and temporal scales, and that model's 
would be greatly augmented by the inclusion of these environmental drivers as regulators of coastal productivity. 
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Introduction 

Climate scenarios suggest that the increased CO2 levels 
in the atmosphere and associated warming will 
increase precipitation in the northern hemisphere 
(Dore, 2005). For northern Scandinavia, precipitation is 
predicted to increase by a maximum of 21 % compared 
to the current average (Meier, 2006). As a result the dis- 
charge of freshwater, and thereby nutrient load to the 
coastal zone, will increase. 

The current view is that increased discharge of nitro- 
gen and phosphorus promotes phytoplankton produc- 
tivity and therefore eutrophication in the coastal zone 
(Larsson et al, 1985; Smith, 2006; Finkel et al, 2010). 
This is based on the fact that most studies implicate 
phosphorus and nitrogen as the major limiting nutri- 
ents of phytoplankton biomass production (Rabalais 
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et al., 2002). Consequently, increased fresh water 
discharge should result in elevated primary production 
and high food web efficiency, and so in greater fish and 
shellfish production (Nixon, 1988). 

However, Howarth et al. (2000), contrastingly reported 
higher primary production during dry than wet years 
in the Hudson estuary, suggesting weaker stratification 
and light penetration to be the main cause. Also, 
organic carbon is typically discharged at higher levels 
simultaneously with nitrogen and phosphorus, poten- 
tially exerting several negative effects on phytoplank- 
ton productivity (Hessen et al., 2010). In accordance 
with this, estuaries worldwide are often reported as 
net-heterotrophic, suggesting that externally supplied 
organic carbon is an important driver of coastal metab- 
olism (Kemp et al, 1997; Sandberg et al, 2004). 
Dissolved organic carbon may support bacterioplank- 
ton biomass production, resulting in increased competi- 
tion for mineral nutrients with phytoplankton 
(Pengerud et al, 1987; Thingstad et al, 2008; Barrera- 
Alba et al, 2009). Increased freshwater discharge may 
further influence stratification of the water column, 
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affecting the vertical distribution of phytoplankton and 
their effective light climate (Cole et al, 1992; Howarth 
et al., 2000; Jager et ah, 2008). Furthermore, an increase 
in humic and suspended substances may reduce light 
climate and, in conjunction with reduced salinity, may 
change the taxonomic composition of resident commu- 
nities (Gasiunaite et ah, 2005; Hessen et al., 2010). 

Marine productivity and transfer efficiency in the 
coastal food web are important to understand in view 
of changing freshwater discharge, as they may influ- 
ence the functionality of the coastal ecosystem. Whether 
larger phytoplankton or small bacterioplankton domi- 
nate the biomass at the food web base will influence the 
food web efficiency, defined as fish production per unit 
of primary production (Nixon, 1988; Rand & Stewart, 
1998; Berglund et al, 2007). If the net effect of increased 
freshwater discharge is higher primary production of 
larger phytoplankton in the coastal zone, this may lead 
to increased production of fish and shellfish (Finkel 
et al., 2010). In contrast, the large influence of riverine 
dissolved and suspended organic matter in the coastal 
zone can both reduce primary production and promote 
a microbial food web structure with poor transfer effi- 
ciency of organic carbon to fish and shellfish (Sandberg 
et al., 2004; Berglund et al., 2007). To properly manage 
coastal resources and mitigate future effects of climate 
change, a better understanding of the net effect of 
increased freshwater discharge on an ecosystem scale is 
therefore required. 

Taken together, it is difficult to predict the net out- 
come of increased river discharge for coastal productiv- 
ity and the balance between auto- and heterotrophic 
processes, or to simulate it in controlled laboratory 
experiments. In this study, we therefore investigated 
the influence of a relatively moderate, but climatologi- 
cally relevant, change in freshwater discharge on the 
ratio of bacterial to phytoplankton biomass production 
in full-scale coastal ecosystems over several years. This 
was used as a proxy for the relative importance of 
microbial production, and thereby losses from fish pro- 
duction and sedimentation. 

Three different coastal sea areas were included, pro- 
viding varying riverine load, morphology, productivity 
level and hydrography. In addition, long-term effects 



were examined by using a 13 year data set. The advan- 
tage of this strategy was that it addressed many of the 
composite effects of increased freshwater discharge on 
hydrography, chemistry, and biology in the coastal zone. 
Long-term field data sets also encompass effects on lar- 
ger and more relevant time scales than can be achieved 
in controlled experimental systems, an important aspect 
in assessing climate-driven environmental change. 

Materials and methods 

Sampling stations 

The study was performed between 1994 and 2006 in three dif- 
ferent basins within the northern Baltic Sea that have varying 
freshwater load, size, and maximum depth (Table 1; Fig. 1). 
Bothnian Bay is divided from the Bothnian Sea by the 20 m 
deep northern Quark (Wulff et al., 1990). The Ore estuary is 
relatively open towards the off-shore Bothnian Sea and lacks a 
marked sill. The catchment area consists mainly of coniferous 
forests and mires, and a smaller proportion of mountainous 
areas with tundra vegetation (Sweden) and arable land 
(Finland). 

Between 13 and 26 samplings per year and basin were per- 
formed from a research vessel, allowing production rate mea- 
surements and sample preparation typically to start within 
1 h from sampling. From 1994 to 1999, one station per basin 
was visited, while sampling was extended to include two rep- 
licate stations in each basin (about 70 km apart) after the 
year 2000. Samples were taken with 7 dm 3 Niskin® bottles 
(General Oceanics Inc., Miami, FL, USA) mounted on a rosette 
sampler. Sampling stations in the Bothnian Bay were A5 
(65°10.00N; 23°14.00E) and A13 (64°42.50N; 22°4.00E); in the 
Bothnian Sea, CI (62°35.20N; 19°58.41E), C3 (62°39.17N; 18° 
57.14E), and C14 (62°5.99N; 18°32.91E); and in the Ore Estuary 
station, B3 (63°29.98N; 19°49.14E) and B7 (63°31.50N; 19° 
48.49E). Station C3 replaced station CI from the year 2000 an 
onwards: both represented the north-western Bothnian Sea. 
Samplings were reduced to 6 yr _1 at station CI 1998, due to a 
temporary reallocation of funding. 

Phytoplankton biomass production (P p ) 

Until the year 2000, water samples were taken at distinct 
depths and incubated at the sampling depth, attached to 
a drifting buoy (Andersson et al., 1996). C0 2 fixation was 



Table 1 Characteristics of sampling stations and basins are shown. Mean load of riverine total organic carbon (TOCr) 1994-2006 is 
presented 
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Fig. 1 The figure shows the study area of the northern Baltic 
Sea, with sampling stations, basins, and bottom topography. 
The Bothnian Bay and Bothnian Sea are divided by the northern 
Quark. The Ore estuary is shown in the inset map. 

measured at five depths by uptake of H 14 CC>3~ 
(2.1 mCi mmor 1 ; Danish Hydrological Institute®) by phyto- 
plankton during 3 h incubation periods, as described in 
Steeman-Nielsen (1952). Depths used were 1, 2, 5, 10, and 
20 m. Daily C0 2 fixation was calculated by multiplying with the 
quotient of integrated diel light irradiance and irradiance dur- 
ing the incubation. Light irradiance in water was measured 
with a Li-Cor LI-193SA spherical sensor (Li-Cor Corporate 
Offices, Lincoln, NE, USA) and surface irradiance with a 
Li-Cor LI-190SA planar sensor (400-700 nm). Trapezoidal 
integration was used to calculate integrated values. From 
2001, samples were collected with a 10 m hose and incubation 
was undertaken onboard the research ship in an ICES incuba- 
tor with differently shaded incubation bottles (Riegman et ah, 
1990). Photosynthetic maximum (P max ) and maximum light 
utilization coefficient (otp : i) was calculated by fitting data to a 
tangential equation according to lassby & Piatt (1976). Daily 
CO2 fixation was calculated from the parameters P max , ot P: i, 
surface reflection of light (empirically derived 0.82) and 
hourly mean irradiance at given depths from light attenuation 
and irradiance below the sea surface. Values down to 20 m 
were modelled using calculated light irradiance at each depth 
and assuming estimated coefficients above 10 m to be applica- 
ble. Light attenuation was determined by fitting In trans- 
formed data to Beer's law, taking deviating surface values, 
ship shadowing, and nonlinear data into account. Trapezoidal 
integration was used to calculate integrated values. Phyto- 
plankton specific growth rate, r p hy to , was calculated as the 
ratio of P p and biomass of phytoplankton. This was measured 
by the Uthermol technique in Lugol preserved samples, using 
a sedimentation chamber according to Andersson et al. (1996). 

A calibration of the different incubation strategies showed 
no significant differences, according to a Model II regression 



of In transformed data (correction factor = 1.08, 95% CI ± 0.27, 
n - 61); neither was a level shift in the time series observed at 
the time of method shift (Fig. SI). 

Bacterioplankton biomass production (P b ) 

Bacterioplankton biomass production rate was measured by 
uptake of 3 H-Thymidine (85 Ci mmoL 1 ; Amersham®; GE 
Healthcare, Buckinghamshire, UK) in 1.5 mL Eppendorf tubes 
(Smith & Azam, 1992). Thymidine uptake in natural samples 
during short incubations (1 h) in small volumes (1 mL) is spe- 
cific for active heterotrophic bacteria, as shown by microautora- 
diography (Fuhrman & Azam, 1982). Thymidine uptake was 
converted to cells produced by an empirically derived factor of 
1.4 x 10 18 cells [mol ^-Thymidine]" 1 (±SE 0.1 x 10 18 ,n = 73), 
calibrating the thymidine uptake to growth of the natural bac- 
terial community (Wikner & Hagstrom, 1999). The carbon den- 
sity per cell was determined by image analysis of acridine 
orange-stained cells and the published volume-to-carbon den- 
sity functions (Norland, 1993; Blackburn et al., 1998). The ratio 
of bacterio- to phytoplankton biomass production rate was 
defined as the Pb:P p ratio. 

Precipitation 

Data on precipitation was provided from the Swedish Meteo- 
rological and Hydrological Institute's regular monitoring pro- 
gram (SMHI, 2009). The station in Haparanda (lat. 65.84, long. 
24.11) was representing the Bothnian Bay and Sundsvall 
airport (lat. 62.52, long. 17.44) represented the Bothnian Sea. 

Riverine flow rate 

Data on river flow was provided from the Swedish Meteoro- 
logical and Hydrological Institute's regular monitoring pro- 
gram (SMHI, 2009). The hydrological database is generally 
based on daily (24 h) data of water flow (m 3 s _1 ). In the stud- 
ied sea area 157 stations, covering 30 major rivers in the drain- 
age area, contributed to the drainage basin aggregated flow 
values. Water flow was measured by water-level gauges 
calibrated by current measurements in river cross-sections 
(Anonymous, 1979; Bergstrom & Carlsson, 1994). For smaller 
watercourses, artificial measuring weirs were used, and dis- 
charge calculated by hydraulic formulas. Monthly means have 
a reported error of ±5%. 

Discharged total organic carbon, total nitrogen, and total 
phosphorus 

Total organic carbon (TOC r ), nitrogen (TN r ), and phosphorus 
(TP r ) in river water were measured by the Swedish University 
of Agricultural Sciences, using standardized protocols within 
the Swedish monitoring program for river run-off (Demandt, 
2009). Within this monitoring program, 25 major rivers in the 
drainage basins were measured at a single sampling station 
per river, per month. TOC r was analysed in HCL-acidified 
samples by high-temperature catalytic oxidation with a 
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Shimadzu TOC r 5050® instrument (Peltzer et al, 1996). Total 
nitrogen and phosphorus were determined by oxidative 
digestion with peroxodisulfate and spectrophotometric analy- 
sis in a Technicon Autoanalyzer (Bran & Luebbe, Solna, 
Sweden), according to the International Organization for 
Standardization (ISO) 11905-1 and ISO 6878, respectively. 

Load of riverine organic carbon to the sea 

The discharge of total organic carbon (TOC r ) was normal- 
ized to the sea basin area to harmonize with biomass pro- 
duction variables in units per area (Table 1). The TOC r was 
assumed to be uniformly distributed in each basin over the 
course of a year. For the Ore estuary, 10% of the TOC r was 
assumed to be metabolized in the estuary, as 90% has been 
shown to be exported to the off-shore areas of the Bothnian 
Sea (Forsgren & Jansson, 1992). Other sources of TOC r to the 
marine environment were assumed to be negligible 
(Sandberg et al, 2004). 

Data aggregation 

Values from distinct depths were integrated by the trapezoi- 
dal method over the average water column in each basin 
(Table 1). Carbon dioxide fixation values were integrated to 
20 m. Where measurements were not available to maximum 
depths, this value was extrapolated from closest relevant 
depths and stations to minimize integration differences. 
Yearly productivity values were calculated from daily values 
by using trapezoidal integrations. 

Statistics 

Data was tested for normal distribution by the Shapiro-Wilks 
test, and for homogeneity of variances according to the Levene 
statistic, using SPSS® statistical software (v. 18; IBM Corpora- 
tion, Armonk, NY, USA). Values were transformed when data 
differed from normal distribution. No significant autocorrela- 
tion was found for discharge of nutrients and freshwater val- 
ues, according to analysis of the autocorrelation function 
(ACF) or partial autocorrelation function (PACF). An a value 
for Type I errors of 0.05 was applied. Tests were two-tailed in 
all cases, where appropriate. Exponential model coefficients 
(cf. Eqn 2) and degree of determination was estimated by the 
curve estimation regression module in SPSS (v. 18). Model II 
linear regression was performed with a reduced major axis 
loss function and bootstrap estimate of coefficient standard 
error. Pearson or Spearman correlations were used, depending 
on the distribution of data. Choice of linear or nonlinear mod- 
els was based on maximizing R 2 and lack of dependence of 
residuals on the predicted values. 

A cubic curve fit was applied using the SPSS® Curve esti- 
mation module to verify that the interannual dynamics with 
time (f) were significant, according to the model 

y = b 0 + bit + b 2 t 2 + b 3 t 3 (1) 

Data sets which lagged 1 and 2 years were analyzed. 
The residuals were tested for autocorrelation in the SPSS® 



Forecasting module and random distribution by graphic plots. 
Where significant autocorrelation was observed, data were 
also applied to an autoregressive integrated model (ARI) with 
appropriate lag and difference of two. Degree of explanation, 
(R 2 ), root mean squared error (RMSE) and other model quality 
measures in the SPSS® Fit syntax were used to assess the best 
model. Cubic curve models showed the best fit to observed 
data, except for Pb:P p vs. time in the Bothnian Sea, where an 
ARI (4,2) model was used. 

A combined standard uncertainty equation was applied to 
estimate the relative contribution of bacterial and phytoplank- 
ton community biomass production, respectively, to the 
change in Pb:P p ratio, as outlined in the published guide to 
the expression of uncertainty in measurement (GUM, 2008, 
Eq. 10). 

Data storage 

All marine data are archived in the regional database dBotnia 
at Umea Marine Sciences Centre (http://www.umf.umu.se) 
and the Swedish national marine database at the Swedish 
Meteorological and Hydrological Institute (SMHI). 

Results 

The annual precipitation showed an elevated level 
between 1998 and 2001 over the drainage and sea area 
for both the Bothnian Bay and Bothnian Sea (Fig. S2). 
The value for year 2000 in the latter area (846.6 mm) 
was the highest precipitation recorded since 1930 
(SMHI, 2009). A natural pulse of freshwater discharge 
consequently occurred during the years 1998-2001 in 
northern Scandinavia, when the discharge was higher 
than in adjacent years (Fig. 2; Table 2). The freshwater 
input was 1.3-1.8 times higher during this high-flow 
period than during 1994-1997 (ANOVA Bonferroni post 
hoc test, P < 0.033 in all cases, n = 4). For the Bothnian 
Sea in the northern Baltic, the flow rates for three out of 
these 4 years were close to the highest recorded in a 
40 year time series, with a peak value 21 % higher than 
the long-term average (91 km 3 yr _1 ) (Bergstrom & 
Carlsson, 1994). 

Total organic carbon (TOC r ) and total nitrogen (TN r ) 
load showed a strong correlation with flow rate, 
explaining the major part of their interannual variation 
(linear regression, R 2 > 0.92, P < 0.001 for all basins, 
n = 13) (On-line data base Swedish University of Agri- 
cultural Sciences). The flow rate also explained a large 
but somewhat smaller part of the variation in total 
phosphorus (TP, R 2 > 0.40, P < 0.021). The variation in 
flow rate in Fig. 2 therefore also indicates the major 
variation in load of TOC r , TN r/ and to some extent TP r 
to the marine environment. 

The load of TN r and TOC r was significantly higher in 
all basins during 1998-2001 than in the preceding 
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Fig. 2 The ratio between bacterioplankton biomass production 
(Pb) and phytoplankton biomass production (P p ) vs. time (filled 
triangles). Non-linear model fits (dashed line) suggested that 
the trends were significant in all cases (R 2 > 0.40, P < 0.037). An 
autoregressive, integrated model was used in the Bothnian Sea. 
River freshwater flow to the basin is also shown (open squares). 
Error bars show the summarized standard deviation of annual 
values based on 5% error of monthly means. The total number 
of samples per year was 72-96 yr . 



period (P > 0.018, ANOVA with Bonferroni's post hoc 
test, Table 2). For TP r , only the P- value for the Bothnian 
Sea (0.097) approached the applied significance level. 
Nonetheless, the mean value of P load during 1998- 
2001 was 7.7%, 21%, and 41% greater than the preced- 
ing period, in the Bothnian Bay, Ore estuary, and 
Bothnian Sea, respectively. The result showed that 
interannual variability in P load was too great to be able 
to detect statistically this magnitude of difference. No 
significant change in the C:N ratio between the periods 
could be shown. However, the C:P ratio of discharged 
material increased during the studied period by about 
50% in the Bothnian Bay and Bothnian Sea (Bonferroni 
test, P = 0.005, n = 4), while the 26% increase in the Ore 
estuary could not be shown to be significant (n = 4). 

We demonstrated a significant shift in the ratio 
between bacterioplankton biomass production and 
phytoplankton biomass production (Pi,:Pp, i.e., 'trophic 
balance') associated with the freshwater pulse (Fig. 2; 
Table SI). A cubic curve fit verified a statistically signif- 
icant peak after year 2000 in the Bothnian Bay 
(R 2 = 0.60, P = 0.034) and the Ore estuary (R 2 = 0.67, 
P = 0.014), while an autoregressive ARI (1, 1) model 
was significant in the Bothnian Sea (R z = 0.40, 



Table 2 Statistics of freshwater and nutrient discharge to the sea basins. Data is aggregated in symmetrical time periods to allow 
statistical tests. Molar ratios of the nutrients are shown. The probability of a Type I error for difference between periods (given by 
the letters) is shown (p). A univariate ANOVA test with Bonferroni's post hoc test for multiple comparisons for equal variances was 
used 



Period mean 



1994-1997 (a) 



1998-2001 (b) 



2002-2006 (c) 



a vs. b 



b vs. c 



Bothnian Bay 

River flow (km 3 yr -1 ) 52 

Total C (10 6 kg yr" 1 ) 258 

Total N (10 6 kg yr- 1 ) 15 

Total P (10 6 kg yr~ J ) 1.0 

CP 678 

C:N 19 

Ore estuary 

River flow (km 3 yr -1 ) 0.85 

Total C (10 6 kg yr _1 ) 8.8 

Total N (10 6 kg yr- 1 ) 0.35 

Total P (10 6 kg yr _1 ) 0.028 

CP 1072 

C:N 29 

Bothnian Sea 

River flow (km 3 yr -1 ) 68 

Total C (10 6 kg yr _1 ) 375 

Total N (10 6 kg yr- 1 ) 22 

Total P (10 6 kg vr -1 ) 0.92 

CP 1085 

ON 20 



70 
440 
24 
1.1 
1020 
22 

1.5 

18 
0.61 
0.034 
1354 

33 

99 
663 
35 
1.3 
1338 
22 



53 
283 

17 
0.74 
978 

19 

0.98 
11 

0.44 
0.022 
1349 
29 

53 
424 
24 
0.70 
1571 
20 



0.033 
0.019 
0.009 
1 

0.005 
0.405 

0.011 
0.022 
0.018 
1 

0.783 
0.263 

0.033 

0.008 

0.007 

0.097 

0.33 

0.36 



0.032 

0.033 

0.24 

0.156 

1 

0.204 

0.03 

0.081 

0.114 

0.661 

1 

0.189 

0.032 

0.017 

0.014 

0.006 

0.36 

0.55 



1 
1 
1 

0.334 
0.008 
1 

1 

1.000 
0.688 
1 

0.738 
1.000 

1.00 

1.00 

1.00 

0.45 

0.016 

0.93 
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P = 0.037). The first derivative at value zero suggested 
the maximum was at year 2002 at all sites (graphic 
determination for the ARI model). A lag of 2 years 
between river discharge and response in Pb:P p ratio 
was indicated in the time series (Fig. 2). Comparing 
1994-1997 to 1998-2001, the relative importance of 
microbial heterotrophic productivity increased in all 
studied sites by a factor of 1.3-2.6. Five years after the 
end of the high flow period (2004-2006), the P b :P p ratio 
was still skewed towards heterotrophy in the Bothnian 
Sea (Table SI, t-test, t = -2.82, df = 4, P = 0.048). In the 
Bothnian Bay and Ore estuary basins, the Pb:P p ratio 
was not significantly different during 2004-2006 com- 
pared to the level before the flow peak. 

Maximum and minimum values of Pb:P p had a ratio 
of between five and seven depending on the site, show- 
ing that variation in trophic balance in the coastal eco- 
system may be substantial (Table SI). However, the 
interannual variability in the Pb:P p ratio appeared stea- 
dier before the peak flow, independent of site 
(±CV = 25%, 1994-1997). Consequently, even a rela- 
tively modest increase in river flow, amounting to 32- 
76% of the 1994-1997 rate (albeit +21% of the long term 
average), had a major effect on the coastal Pb:P p ratio in 
our case study. 

The importance of microbial heterotrophic processes 
(Pb:P p ) increased exponentially with increased dis- 
charge of riverine TOC to the coastal zone (TOC r ), 
when lagged +2 years according to the model 

P b :P p = ax e bxTOC ' (2) 

This was valid for all sites, based on statistically sig- 
nificant coefficients from an exponential regression 
(Fig. S3, Table 3). TOC r explained a significant part of 
the variation (i.e., R 2 ) in Pb:P p at all sites, with the 



Table 3 Statistics of exponential curve fit regression for the 
relationship between P b :P p and specific TOCr discharge 
(lag = year + 2). The Y-intercept (a), exponential coefficient 
(b) and their standard error (SE), probability of type I error 
(P), and number of annual values analyzed (n) are shown. 
One missing P b :P p value in each off-shore basin was interpo- 
lated 

(-) ([g CP 1 m 2 yr) 

Basin a SE (a) P (a) b SE (b) P (b) n 

Bothnian 0.56 0.23 0.035 0.093 0.041 0.049 11 
Bay 

Ore 0.36 0.12 0.015 0.0050 0.0016 0.013 11 
estuary 

Bothnian 0.26 0.17 0.15 0.18 0.081 0.053 11 
Sea 



highest value at the coastal site. The y-intercepts were 
significant in the Bothnian Bay and Ore estuary. Only 
the b coefficient in the Ore estuary was different from 
the off-shore basins, based on 2 x SE, a consequence 
of the higher TOC r to this coastal site. A model II linear 
regression also showed significantly positive slope 
coefficients (data not shown). 

Most of the change in Pb:P p was due to a decrease in 
phytoplankton biomass production, rather than an 
increase in bacterioplankton biomass production 
(Fig. SI). A variance component analysis (standard 
uncertainty calculation) suggested that, on average, 
71% of the variation in the Pb:P p ratio was due to the 
decrease in phytoplankton biomass production and 
29% was due to increase in bacterioplankton biomass 
production (Table S2). The difference between mini- 
mum and maximum phytoplankton annual production 
was also unexpectedly large, threefold to sixfold 
depending on sea basin (Fig. SI). However, where 
more freshwater discharge promoted increased input 
of TN r and TP r , and thus potentially increased phyto- 
plankton production, our data showed that it instead 
decreased (Fig. SI; Table 2). 

To test for a potential effect on phytoplankton due 
the change in growth conditions, r ph y to was analyzed 
with specific TOCr discharge by a model II regression 
and Pearson correlation (Fig. 3). This showed a signifi- 
cant negative relationship for both the Ore estuary 
(P < 0.05, R 2 = 0.46) and the Bothnian Sea (P < 0.05, 
R 2 = 0.32). The Bothnian Bay did, however, not show 
a significant correlation, while a negative slope 
coefficient was demonstrated (P < 0.05). 

Discussion 

Increased precipitation around the millennium shift 
caused a marked increase in freshwater discharge dur- 
ing a few years, and also direct input of freshwater and 
nutrients to the sea surface by rain (Figs 3 and 4). By 
using this extreme peak flow of river water in a 40 year 
time series, in conjunction with long-term ecological 
marine time series data, we demonstrated a marked 
shift from autotrophy (Pb:P p < 1) to heterotrophy (P b : 
P p > 1) in three coastal sites (Fig. 2; Table SI). At all 
sites, the Pb:P p ratio was positively correlated with the 
load of TOC r (Fig. S3; Table 3). This result showed that 
the net effect of increased freshwater flow to the marine 
environment hampered phytoplankton biomass 
production, while bacterial biomass production was 
maintained. This observation was unexpected consider- 
ing that increased river flow was accompanied by an 
increase in TN r and TP r , implying increased phyto- 
plankton biomass production on mineralized nutrients. 
However, the observed effect may be explained 
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Fig. 3 Phytoplankton specific growth rate (r p hy to ) is shown as a 
function of the specific riverine TOC discharge. The line of natu- 
ral logarithm transformed values was derived using a Model II 
linear regression with a major axis loss function. The slope coef- 
ficient for the Bothnian Bay had a P value > 0.05. 

by including influence of freshwater, TOC r , and 
bacterioplankton growth in production models for the 
coastal ecosystem, as outlined below. 

The discharge of riverine TOC was a major driver of 
the Pt,:P p ratio at all sites, assuming some contribution 
to the remaining variation from random measurement 
variance (R 2 - values, Fig. S3). The slope coefficients 
were not significantly different between the off-shore 



basins, suggesting a similar response on trophic balance 
from increased TOCr discharge. The markedly lower 
slope coefficient for the coastal Ore estuary reflects the 
markedly higher specific load of riverine TOC at this 
coastal site. However, due to the correlation between 
TOC r and freshwater discharge {R 2 > 0.92), either factor 
may still contribute to explain the observed change in 
trophic balance. The degree of explanation was slightly 
lower in the off-shore basins when using freshwater 
discharge as the independent variable, and the slope 
coefficient for the Bothnian Sea was not significantly 
different from zero. This indicates a specific importance 
of discharged TOCr for the observed effects. Other vari- 
ables investigated to explain the variability in the Pt,:P p 
ratio, such as the C:P ratio in river discharge, marine 
nutrient concentrations or global photosynthetically 
active irradiation, did not show significant correlation 
in any basin (data not shown). The only exception was 
marine inorganic nitrogen, which showed a positive 
relationship with Pb,:P p in the Bothnian Sea (r = 0.80, 
P = 0.02). 

Our results suggest that the suppression of phyto- 
plankton biomass production explained a greater pro- 
portion (71%) of the change in Pb:P p ratio than an 
increase in bacterial biomass production (Table S2). 
Therefore, an explanation of the change in Pt,:P p ratio 
should be primarily sought in the control of phyto- 
plankton biomass production. One contributing reason 
for the lack of a positive response in phytoplankton bio- 
mass production was that increased freshwater flow 
mainly increased the TOC r and TN r load, while an 
increase in TP r was moderate and not statistically 
significant (Fig. SI; Table 2). Thus, the increase in the 
major limiting nutrient, phosphorus, was lower relative 
to other nutrients, at elevated flows of freshwater. 
Phosphorus is considered the limiting nutrient in most 
of this study area, while a shift to nitrogen limitation 
occurs in the southern Bothnian Sea (Zweifel et al., 
1993; Andersson et al., 1996). Total nutrient and carbon 
pools were used, as these also encompass organic and 
particulate fractions, assumed to be largely available 
also to osmotrophs by phagotrophy, mineralization, 
and photolysis in the food web during the long time 
scales investigated (Kieber et al., 1990; Sandberg et al., 
2004). 

The lag applied in Table 3 means that the TOCr dis- 
charge at a given year resulted in a full effect showing 
in the trophic balance 2 years later (Fig. 2). Based on 
the primary data, this is caused by the reduced bacte- 
rial biomass production in the first year of high fresh- 
water discharge (year 1998, Fig. SI). The P b then 
successively increased to preflow level or above, also 
elevating the Pt,:P p ratio. A hypothesis to explain the 
observed lag period may be that initially recalcitrant 
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carbon and nutrient substrates (organic and particle 
bound), are gradually made available by exudates 
released during phagotrophy of particles, mineraliza- 
tion, and photolytic cleavage of dissolved organically 
bound nitrogen and phosphorus, when exposed to the 
marine environment (Pengerud et al., 1987; Kieber 
et al, 1990; Thingstad et al, 1993). In addition a change 
in bacterioplankton diversity, as an adaptation to the 
newly available substrate resources, may gradually 
promote a more efficient metabolism of the terrigenous 
substrates, based on observations in culture studies 
(Kisand et al, 2002). The response of the phytoplankton 
specific growth rate (i.e., r phyto ), in contrast, was imme- 
diate, as shown by a significant relationship with TOCr 
discharge, without lag, discussed below (Fig. 3). 

Field data supported the hypothesis that growth limi- 
tation of phytoplankton was present because of a signif- 
icant negative relationship between )' phyto and specific 
TOCr discharge (Fig. 3). The realized r phyto is expected 
to decrease with enhanced competition, and thus, in 
this case, with increased Pt,:P p . The strongest influence 
of TOC r (i.e., highest slope coefficient) occurred in the 
Bothnian Sea, suggesting this basin to be most sensitive 
to elevated freshwater discharge. The negative effect on 
'"phyto could, however, also result from other factors, 
e.g., reduced effective light climate as a consequence of 
discharged coloured and suspended matter, or changed 
stratification (Howarth et al, 2000; Hessen et al, 2010). 
This was based on that the results were almost identical 



when using freshwater discharge as the independent 
variable (data not shown). We can only speculate that 
the lack of a significant correlation in the Bothnian Bay 
was due to the smaller contribution of phytoplankton 
community growth to the observed change in trophic 
balance in this basin (Table S2). 

The effect of freshwater discharge on phytoplankton 
growth without lag suggested an influence of factors in 
the discharge which exerted an immediate action 
(Figs 3 and SI). Discharge of coloured humic substance, 
suspended material and a deeper mixed layer, reducing 
the effective light climate for phytoplankton, are all in 
synergy and expected to act without a time lag on the 
annual scale (Cole et al, 1992; Howarth et al, 2000; 
Hessen et al, 2010). This explanation is also in accor- 
dance with model analyses, suggesting light limitation 
to override nutrient limitation in deeper and better 
mixed water columns (Jager et al, 2010). The simulta- 
neously increased load of carbon to the coastal zone 
provided another, synergistic, explanation for the shift 
in carbon balance, and reduction in phytoplankton bio- 
mass production. Riverine dissolved organic carbon is 
an alternate carbon source of phytoplankton released 
carbon for the heterotrophic bacteria (Zweifel et al, 
1993; Wikner et al, 1999). This may have promoted a 
competition between bacterio- and phytoplankton for 
inorganic nutrients, where the former has superior 
affinity for these substrates (Pengerud et al, 1987; 
Mindl et al, 2005). As bacterioplankton community 
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growth showed a lag of 2 years, competition for limit- 
ing nutrients is, however, not likely to explain the ini- 
tial reduction in phytoplankton growth (Fig. SI). When 
the bacterioplankton community with time became bet- 
ter able to exploit the more riverine-influenced carbon 
source the competition scenario may have contributed 
by adding a synergistic hampering effect on phyto- 
plankton production. Simultaneous discharge of 
elements such as Fe associated with TOCr discharge 
may in addition, by chemical binding, have reduced 
the available pool of P0 4 ~ for phytoplankton (Blomq- 
vist et al, 2004). Consequently, reduced access to min- 
eral nutrients may also explain the hampered 
phytoplankton biomass production. The relative impor- 
tance of changed light conditions for phytoplankton 
and competition of limiting nutrient was however not 
investigated in this study. 

Similar studies of the composite effects of discharged 
carbon, nitrogen, and phosphorus on trophic balance in 
the coastal zone are scarce in the literature, especially 
given the datasets cover large temporal and spatial 
scales. In a long-term study of the Hudson estuary, Ho- 
warth et al. (2000) stressed the importance of climate 
drivers, such as freshwater discharge, on primary pro- 
duction. Similar to our findings, primary production 
was suppressed at higher levels of freshwater discharge 
in their study, primarily explained by weaker stratifica- 
tion and light penetration. Caraco & Cole (2003) also 
reported a tenfold lower ratio of de novo gross primary 
production to allochthonous organic load during a wet 
year than a dry one in the Hudson River, in agreement 
with our conclusions. In an estuary-lagoon system in 
eastern Brazil, increased freshwater flow was found to 
stimulate bacterial production while hampering pri- 
mary production, with a resulting Pb:P p ratio of up to 
20 during the rainy season (Barrera-Alba et al, 2009). 
Our findings were also in line with the Thingstad-para- 
dox' observed in an Arctic mesocosm experiment, in 
which addition of carbon substrate led to a lower accu- 
mulation of carbon in the system, due to hampered 
phytoplankton biomass production (Thingstad et al., 
2008). Similar results are reported from the Hunter 
Estuary, Australia, where addition of glucose reduced 
chlorophyll-s concentration, while increasing bacterial 
abundance and respiration (Hitchcock et al., 2010). In 
shallow unproductive lakes increased load of terrige- 
nous dissolved organic carbon (DOC) is also responsi- 
ble for hampering total lake productivity (Karlsson 
et al., 2009). Thus, our results add support to an emerg- 
ing picture of TOC availability as an important, but 
counterintuitive, controller of aquatic productivity. This 
study contributes by showing the importance of 
riverine TOC for the coastal zone at large temporal and 
spatial scales. It also provides evidence for a general 



relationship between the riverine TOC supply and P b : 
P p ratio of the coastal sea. 

The explanatory model outlined above seems mainly 
associated with estuaries at temporal scales from sea- 
sonal to annual. In short-term studies of freshwater 
plumes entering directly into oceanic waters (i.e., Strait 
of Georgia, Bay of Bengal), a stimulating effect on pri- 
mary production has been reported (Harrison et al., 
1991; Gomes et al., 2000). Effects were, however, 
reported to alternate between negative and positive, 
depending on the temporal scale in one particular envi- 
ronment. Thus, in oceanic environments positive effects 
of freshwater plumes on plankton productivity are pos- 
sible, suggesting a lower negative influence of the 
freshwater and TOC supply on phytoplankton growth 
in these environments. 

We consider the bacterial and phytoplankton bio- 
mass production ratio to have sufficient accuracy when 
estimating the distribution of major pelagic carbon 
flows (i.e., microbial hetero- vs. autotrophic processes). 
Phyto- and bacterioplankton biomass production 
clearly dominates the marine food web carbon flow, 
and are therefore pivotal factors in coastal ecosystem 
function (Kemp et al, 1997; Sandberg et al, 2004). P b :P p 
values have further been shown to give similar results 
to alternative methods using mass balances or oxygen 
models in estimating net ecosystem metabolism, net 
CO2 emission and net heterotrophy (Stigebrandt, 1991; 
Kemp et al, 1997; Sandberg et al, 2004). It is also well 
documented that increase in bacterial biomass produc- 
tion will promote carbon flow through the multi- 
trophic level protozoan dominated food chain, with 
associated respiration losses (Pengerud et al, 1987; 
Wikner & Hagstrom, 1988; Sandberg et al, 2004; Bergl- 
und et al, 2007). Lower fish production and sedimenta- 
tion is thereby expected, compared to the shorter food 
chain associated with the flow of phytoplankton bio- 
mass directly through the mesozooplankton commu- 
nity. Bacterio- and phytoplankton biomass production 
are therefore key processes for the trophic balance and 
food web efficiency of the coastal sea. 

The found relationship between riverine discharge of 
TOC and coastal trophic balance suggests that 
increased precipitation, and the resulting higher load of 
freshwater and terrigenous TOC, will lead to a food 
web dominated by heterotrophic microorganisms 
(Sandberg et al, 2004; Berglund et al, 2007; Fig. 4). As 
production of organic matter by phytoplankton is 
simultaneously suppressed, a lowered supply of food 
to fish and benthic organisms is expected (Nixon, 1988; 
Berglund et al, 2007; Finkel et al, 2010). Our view of 
the effect of precipitation on coastal productivity is thus 
opposite to recent assessments of climate-change effects 
for the Baltic Sea (Neumann, 2010). Neumanns' s model 
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study suggested an increase in phytoplankton biomass 
of 5% and extended growth season at elevated precipi- 
tation based on the classical view of phytoplankton 
growth limited only by phosphorus or nitrogen. Influ- 
ence of freshwater discharge on pelagic stratification, 
light climate, and concentration of riverine TOC should 
therefore be included in current models for coastal 
management and resource utilization (e.g., mitigation 
of eutrophication and allowable fish catch), as well as 
in forecasting the consequences of climate change 
(Fig. 4). Supply of nitrogen and phosphorus alone are 
insufficient to accurately predict coastal productivity. 
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